† Corresponding author. E-mail:
Project supported by the National Natural Science Foundation of China (Grant Nos. 11374217 and 11176020).
Ab initio calculations of lattice constants, lattice stabilities of HgX (X = S, Se, Te) at different electronic temperatures (Te) have been performed within the density functional theory (DFT). We find that the lattice constants of HgX increase and the phonon frequencies reduce as Te increases. Especially the transverse-acoustic (TA) phonon frequencies of HgX gradually become negative with the elevation of the electron temperature. That is to say ultrafast intense laser induces lattice instabilities of HgX and athermal melting appears for the increase of laser intensity. What is more, with the X atom number increasing, the critical electronic temperatures of HgX are decreased in sequence. This result would be helpful for understanding the athermal melting processes for femtosecond laser micromachining.
The mercury chalcogenides HgX (X = S, Se, Te) are attracting increasing attention, because of their technological interest for electro–optical devices and spintronic applications.[1,2] As is known, HgS is a technologically important semiconductor for dichorism,[3] birefringence,[4] photoelectric,[5,6] and acousto–optic properties.[7] Recently, the structural, electronic, and thermal properties of HgX binary compounds have been investigated within the density functional theory (DFT).[8–14] Al Shafaay et al.[8] applied a first principle full potential-linearized augmented plane wave (FP-LAPW) method to study the structural, electronic, and thermal properties of HgX binary compounds as well as HgSxSe1−x and HgSexTe1−x ternary alloys. Svane et al.[9] have shown that HgS is a semiconductor with gap energy equal to 0.61 eV using hybrid quasi-particle selfconsistent GW (QSGW) approximations. Cardona et al.[10] also used the hybrid QSGW approach to investigate the electronic band structure and the phonon dispersion relations of the zinc blende type mercury chalcogenides, which showed that HgS is a semiconductor, and both HgSe and HgTe have inverted band structures with Γ6 lying below Γ8. Besides, Verma et al.[12] calculated the thermal properties of mercury chalcogenides. Secuk et al.[13] investigated the structural optimization, electronic band structure, density of electron states, optical, dynamic and thermodynamic properties of HgTe and HgSe within local density approximation (LDA). Wang et al.[14] reported that all II–VI semiconductors with a cubic zinc blende structure have negative thermal expansion (NTE) behavior at low temperatures within the framework of first-principles.
Material response under intense laser irradiation is of significant importance for various applications, such as laser micromachining, laser induced damage, surface treatment, and surgery.[15] Therefore, the interplay between ultrafast laser and material requires a comprehensive theoretical description. A large number of theoretical and experimental efforts have thus been stimulated to study the ultrafast laser–materials interactions.[16–22] Recent advances have revealed that semiconductors under intense laser irradiation firstly underwent an athermal transition induced by the softening inter-atomic bonds rather than conventional thermal transition.[16–20] Feng et al.[16] have investigated the photo-induced non-thermal phase transition of β-SiO2 by ab initio calculation. To clarify the electronic-excitation effect under intense laser illumination, Recoules et al.[17] employed density functional perturbation theory (DFPT) to compute the evolution of the phonon spectrum at different electronic temperatures for Si, Al, and Au, which provides a comprehensive study of the lattice stability considering the effect of the laser illumination, meanwhile the electronic excitation in silicon can induce lattice instabilities as exhibited by transverse-acoustic (TA) phonon modes with imaginary frequencies. Shen et al.[18] studied the effect of intense laser irradiation on the lattice stability of β-SiC by ab initio calculation within the density functional perturbation theory (DFPT).
In the case of ultrashort laser irradiation, semiconductors produce a large number of strong inhomogeneous heating electrons. The heating electrons in the semiconductors rapidly exchange energy through electron–electron collisions and then instantaneously reach a high temperature (thousands of Kelvins, electronic temperatures) on the tens of femtosecond timescale.[23] Accurate knowledge of the rapidly evolving conditions is of importance.[24] At the same time, the inter-atomic potentials of semiconductors encounter dramatic changes and the lattices destabilize then the solids catch ultrafast structural transitions. A few groups have confirmed that the laser-semiconductor interaction occurs in an athermal transition.[16,17,22] Therefore, accurate evaluation of the electronic excitation is essential not only from a physical point of view but also for device engineering.
Even though there has been remarkable experimental and theoretical work on HgX, there are limited theoretical studies on the effect of ultrafast intense laser irradiation on HgX in the literature. It is very necessary to investigate the lattice stability of HgX considering the electronic excitation effect through ab initio calculations based on DFT. In this work, we aim to employ ab initio pseudopotential method to study the lattice constants, transverse-acoustic phonons of HgX at different electronic temperatures caused by different intensities of ultrafast intense laser irradiation. Finally, the lattice stabilities of HgX are further discussed at higher electronic temperatures.
The binary compounds HgX (X = S, Se, Te) crystallize in the zinc-blende structure and their space groups are F-43m, meanwhile the Hg atom is located at the origin and the X atom is located at (0.25, 0.25, 0.25).
We performed all calculations based on finite-temperature density functional theory (DFT) on the lattice parameters and dynamic properties of HgX with electronic temperatures (Te) increasing from 0 eV to 4.0 eV. For the treatment of exchange and correlation, the generalized gradient approximation (GGA) parameterized by Perdew–Burk–Ernzerhof (PBE)[25–27] is used. The 3s2 3p2, 4s2 4p2, 5s2 5p2, and 5d10 6s2 are retained as valence electrons for S, Se, Te, and Hg, respectively. The projector-augmented wave (PAW)[28] method is employed to describe the electron–ion interactions as implemented in the Vienna ab initio simulation package (VASP).[29–31] The plane-wave kinetic energy cutoff is set to 370 eV for HgS, and 310 eV for HgSe and HgTe.[10] For the sampling of the Brillouin zone, we use 11 × 11 × 11 k-point grids for geometry optimization at different electronic temperatures. All parameters have been tested for convergence. Lattice parameters and internal positions at various electronic temperatures in primitive cells are fully relaxed until residual Hellman–Feynmann (HF) forces are less than 1.0 × 10−4 eV/Å and the energy convergence threshold is set as 1.0 × 10−7 eV/atom. Moreover, a few hundred bands are chosen to make enough electrons occupy the conduction band at high electronic temperatures.
Phonon-dispersion calculations are performed by inter-atomic force constants in the real space[32,33] within density functional perturbation theory (DFPT).[34] In the calculation we use the supercell approach to achieve the real space force constants together with 16 atoms 2 × 2 × 2 supercells, a Γ-centered 5 × 5 × 5 k mesh. The phonon dispersion curves of HgX compounds in the considered structures are calculated by using the PHONOPY code[35] based on the frozen phonon method. Since HgX belongs to an ionic crystal, atomic displacement induces a dipole and inter-atomic force constants are remarkably influenced by dipole–dipole interaction which causes longitudinal optical–transverse optical (LO–TO) splitting at the wavelength vector of
For the zinc-blende HgX, the calculated lattice constants a0 at the ground state, compared with experimental and other theoretical results, are listed in Table
![]() | Table 1. Lattice constants of zinc-blende HgX are obtained by the ground-state calculation compared to other theoretical and experimental studies. . |
Figure
To explore ultrafast intense laser inducing the lattice stabilities, the phonon-dispersion curves at four electronic temperatures are calculated, as shown in Fig.
Figure
In addition, we present TA phonon frequencies at high symmetry L point as a function of Te in Fig.
For zinc-blende prototype II–VI semiconductors, the TO–LO phonon frequencies at the Γ point can be expressed by the following Lyddane–Sachs–Teller (LST) relation:[43] ωLO/ωTO = [ε(0)/ε(∞)]1/2, where ε(0) and ε(∞) represent the static and high-frequency dielectric constants, respectively. As ionic crystals, the effective charge amount of positive and negative ions can affect the strength of the polarization electric field. The more effective charge amounts of ions attribute to the bigger divisions between ωLO and ωTO.[43] Consequently, the dipoles induce the electric field and influence the LO–TO splitting degree (δω) which can assess the ionic strength of ionic crystals at the Γ point.[19] Seen from Table
![]() | Table 2. Experimental and calculated phonon frequencies of HgX for the high symmetry point Γ in the Brillouin zone, in units of cm−1. . |
The LO–TO splitting degrees (δω) of HgX (X = S, Se, Te) at the Γ point are a function of the Te as shown in Fig.
In conclusion, we have studied the electronic excitation effects of ultrafast intense laser exposure on the zinc-blende mercury chalcogenides HgX (X = S, Se, Te) at different electronic temperatures based on DFT. It is indicated that the lattice constants of HgX increase and all the phonon modes soften as the Te is increased. The crystal lattices of HgX become more instable through the more negative TA phonon frequency due to the thermal electronic excitations. When the electronic temperature increases to 1.45 eV for HgS, 1.31 eV for HgSe, 1.20 eV for HgTe, the entire TA frequencies become negative, which gives signatures of athermal phase transitions. With the X atom number increasing, the critical electronic temperatures successively decrease. Moreover, as for semimetal HgSe and HgTe, the LO–TO splitting degrees (δω) are different from that of semiconductor HgS. This result can contribute to a good understanding for athermal melting processes of the femtosecond laser micromachining devices and materials.
1 | |
2 | |
3 | |
4 | |
5 | |
6 | |
7 | |
8 | |
9 | |
10 | |
11 | |
12 | |
13 | |
14 | |
15 | |
16 | |
17 | |
18 | |
19 | |
20 | |
21 | |
22 | |
23 | |
24 | |
25 | |
26 | |
27 | |
28 | |
29 | |
30 | |
31 | |
32 | |
33 | |
34 | |
35 | |
36 | |
37 | |
38 | |
39 | |
40 | |
41 | |
42 | |
43 | |
44 | |
45 |